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Flowering time in plants is a tightly regulated process. In barley (Hordeum vulgare L), 
HvFTl, ortholog of FLOWERING LOCUS T, is the main integrator of the photoperiod and 
vernalization signals leading to the transition from vegetative to reproductive state of the 
plant. This gene presents sequence polymorphisms affecting flowering time in the first 
intron and in the promoter. Recently, copy number variation (CNV) has been described 
for this gene. An allele with more than one copy was linked to higher gene expression, 
earlier flowering, and an overriding effect of the vernalization mechanism. This study aims 
at (1) surveying the distribution of HvFTl polymorphisms across barley germplasm and (2) 
assessing gene expression and phenotypic effects of HvFTl alleles. We analyzed HvFTl 
CNV in 109 winter, spring, and facultative barley lines. There was more than one copy 
of the gene (2-5) only in spring or facultative barleys without a functional vernalization 
VrnH2 allele. CNV was investigated in several regions inside and around HvFTl. Two 
models of the gene were found: one with the same number of promoters and transcribed 
regions, and another with one promoter and variable number of transcribed regions. This 
last model was found in Nordic barleys only. Analysis of HvFTl expression showed that 
association between known polymorphisms at the HvFTl locus and the expression of 
the gene was highly dependent on the genetic background. Under long day conditions 
the earliest flowering lines carried a sensitive PpdHl allele. Among spring cultivars with 
different number of copies, no clear relation was found between CNV, gene expression 
and flowering time. This was confirmed in a set of doubled haploid lines of a population 
segregating for HvFTl CNV. Earlier flowering in the presence of several copies of HvFTl 
was only seen in cultivar Tammi, which carries one promoter, suggesting a relation of gene 
structure with its regulation. HvCEN also affected to a large extent flowering time. 
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INTRODUCTION 

In temperate cereals, like barley (Hordeum vulgare L.) and wheat 
(Triticum aestivum L.), flowering is regulated by the integra- 
tion of two seasonal signals (Laurie, 2009): photoperiod (day 
length) and vernalization (prolonged exposures to low temper- 
atures). Flowering time is also closely linked with agronomic 
performance. Plants must flower at the appropriate time of the 
year, when conditions are most favorable for pollination, seed 
development and high grain yield. 

The responses to day length and temperature serve to classify 
barley varieties according to their adaptation pattern. Based on 
the response to day length, varieties are divided into photoperiod- 
sensitive (long days accelerate flowering) or -insensitive (plants 
flower almost independently of the day length). Based on the 
response to vernalization, barley varieties are classified as win- 
ter (vernalization is required for timely flowering) or spring 
(flowering irrespective of vernalization), although the presence 
of an allelic series at VrnHl produces intermediate genotypes 
(Hemming et al., 2009). Usually, winter varieties are sown in 
autumn, spring varieties in winter and spring, and there is a third 
category known as facultative varieties, that can be sown anytime. 



Several major genes are the main responsible for the responses to 
photoperiod and vernalization, which are described next. 

Allelic differences in the photoperiod genes PpdHl and PpdH2 
are associated with natural variation in the response to day length. 
PpdHl (a member of the Pseudo Response Regulator family) is 
part of the circadian clock of the plant and promotes flowering 
under long days (Turner et al., 2005). Recessive mutations in the 
PpdHl gene result in delayed flowering under long days (Turner 
et al., 2005; Hemming et al., 2008). The PpdHl gene acts in par- 
allel to HvCOl (Campoli et al, 2012), which is one of the barley 
homologs of the Arabidopsis (Arabidopsis thaliana) CONSTANS 
(CO) gene. Overexpression of HvCOl results in the up-regulation 
of HvFTl (the ortholog in barley of the Arabidopsis FLOWERING 
LOCUS T, or FT) and the acceleration of flowering (Campoli 
et al., 2012). PpdH2 has been long acknowledged as the responsi- 
ble for acceleration of flowering in response to short photoperiod, 
although its role is being re-defined (Casao et al, 201 1). It is a par- 
alog of HvFTl, (alternatively named HvFT3 by Faure et al, 2007 
and Kikuchi et al., 2009). Its effect on flowering is not as strong 
as HvFTl and it seems to be restricted to winter genotypes under 
short days or long days without vernalization (Casao et al, 2011). 
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Another paralog of HvFTl with a large effect on time to flowering, 
particularly at Mediterranean latitudes (Boyd et al, 2003; Cuesta- 
Marcos et al., 2008), is HvCEN, and its two main haplotypes are 
differentially distributed over spring and winter barley varieties 
(Comadran et al., 2012). 

Natural variation in barley vernalization requirement is mainly 
found in the vernalization loci VrnHl (Trevaskis et al., 2003; Yan 
et al, 2003), VrnH2 (Yan et al, 2004), and VrnH3 (Yan et al., 
2006). The VrnHl gene is closely related to the Arabidopsis gene 
APETALA1 (API), responsible for the transition from the vege- 
tative to the reproductive stage (Trevaskis et al, 2007). Different 
alleles of VrnHl have been identified, with insertions or deletions 
in the first intron of the gene (von Zitzewitz et al., 2005; Cockram 
et al., 2007; Hemming et al, 2009), affecting the length of the 
optimum vernalization period. Alleles lacking large sections of 
the ~11 kb intron are more active and are associated with ear- 
lier flowering without vernalization, whereas alleles lacking small 
segments are associated with only a moderate increase in VrnHl 
activity and weaker promotion of flowering (Szucs et al, 2007; 
Hemming et al., 2009). VrnH2 includes three closely related genes 
designated as HvZCCTa-c, which are characterized by a putative 
zinc finger and a CCT-domain. VrnH2 is considered to play the 
role of repressor of flowering and it has been shown that dele- 
tions of all the three HvZCCT genes result in spring growth habit 
(Karsai et al., 2005; Trevaskis et al., 2006). Finally, VrnH3 was 
shown to correspond to HvFTl, the ortholog of the Arabidopsis 
FT gene (Yan et al, 2006; Faure et al., 2007). 

FT is considered as the main flowering integrator of the pho- 
toperiod and vernalization pathways in both monocot and dicot 
species (Turck et al., 2008). In barley, the most accepted hypoth- 
esis for the regulation of HvFTl establishes that, during the fall, 
when temperate cereals germinate, HvFTl is repressed by VrnH2 
(Yan et al., 2006; Hemming et al, 2008; Distelfeld et al, 2009). 
During winter, vernalization up-regulates VrnHl (Trevaskis et al., 
2006; Oliver et al., 2009), which results in the repression of VrnH2 
in the leaves and, consequently, the activation of HvFTl transcrip- 
tion in the spring (Loukoianov et al., 2005; Trevaskis et al, 2006; 
Hemming et al., 2008; Chen and Dubcovsky, 2012). The precise 
molecular mechanisms of action of these genes are still the object 
of numerous studies in barley and other cereals. FT transcription 
is induced in the leaves and it has been demonstrated in differ- 
ent species that the encoded protein travels through the phloem 
to the stem apical meristem, where it plays a central role in trig- 
gering flowering (Corbesier et al, 2007; Tamaki et al., 2007). In 
Arabidopsis, FT interacts with the bZIP transcription factor FD 
and up-regulates the expression of the meristem identity gene 
API at the shoot apex (Abe et al, 2005; Wigge et al., 2005). A simi- 
lar interaction has been reported in wheat, where the orthologous 
FT protein interacts with an FD-like protein and has the ability to 
bind in vitro the promoter of VrnHl, the wheat homolog of API 
(Li and Dubcovsky, 2008). 

Ample natural variation in the HvFTl gene has been found, 
with polymorphisms reported in the promoter and in the first 
intron. This variation has been linked to differences in flow- 
ering phenotypes in a number of studies (Yan et al., 2006; 
Hemming et al., 2008; Casas et al., 2011). It seems clear now 
that, in Arabidopsis, the FT promoter and first intron contain 



cis-regulatory sites that are important for its transcriptional reg- 
ulation (Tiwari et al., 2010). However, the FT1 regulatory regions 
of barley and wheat are not as well characterized. Yan et al. 
(2006) found an association between growth habit and muta- 
tions in the first intron, but the sequencing of the HvFTl alleles 
from populations previously used to map QTL for flowering time 
(Hemming et al., 2008) failed to reveal any significant associa- 
tion between the two linked SNPs in intron 1 and flowering time. 
Further results reported in other surveys of HvFTl allelic varia- 
tion (Cuesta-Marcos et al., 2010; Casas et al., 2011) were also in 
disagreement with Yan et al. (2006) regarding the direction of the 
effect assigned to the functional polymorphism in the first intron. 
Yan et al. (2006) also identified two promoter haplotypes, charac- 
terized by seven linked SNPs and two indels in the first 550 bp 
upstream of the start codon. Using primers specific to differen- 
tiate those indels, Casas et al. (2011) analyzed natural variation 
for the promoter haplotypes (135-146 vs. 139-142 bp) and the 
intron 1 haplotypes (AG/TC) in a collection of barley landraces. 
In that study the intron AG haplotype was clearly associated with 
later flowering than the TC haplotype. The results for the pro- 
moter haplotypes hinted at a role of these polymorphisms on 
flowering time, but of lesser magnitude than intron polymor- 
phism. The combination of the 135-146 promoter with the TC 
intron was associated with earliest flowering (Casas et al., 2011; 
Ponce-Molina et al, 2012). Further evidence from other popula- 
tions (Nitcher et al., 2013) confirmed the description of promoter 
haplotypes as "early" (135-146) and "late" (139-142). Another 
SNP in the promoter of the HvFTl gene, upstream of the studied 
region was also suggested to have an additional role on flowering 
time (Cuesta-Marcos et al., 2010; Casas et al, 2011). 

The scope of HvFTl polymorphism has been recently widened 
even further by including copy number variation (CNV), first 
described by Nitcher et al. (2013). Recently, this type of polymor- 
phism has been proposed as a key contributor to intra-species 
genetic variation, along with SNPs and indel polymorphisms. 
Some data suggest that CNV mainly affects the members of 
large families of functionally redundant genes and that the effects 
of individual CNV events on phenotype are usually modest 
(Zmienko et al., 2014). Nevertheless, there are many cases in 
which CNVs for specific genes have been linked to important 
traits such as flowering time and plant height and resistance 
(Zmienko et al., 2014). Regarding HvFTl, it has been reported 
recently that a genotype with high gene copy number (BGS213, 
derived from cultivar Tammi) was responsible for early flowering 
and an epistatic override of winter growth habit, caused by the 
combination of vrnHl and VrnH2 alleles (Nitcher et al., 2013). 
This study pursues to carry out a comprehensive survey of HvFTl 
polymorphisms, including CNV, in barley accessions of different 
origins and germplasm groups. Also, we aim to provide further 
new information on gene expression and phenotypic effects of 
contrasting genotypes at HvFTl. 

MATERIALS AND METHODS 
PLANT MATERIALS 
HvFTl polymorphism 

A set of 109 genotypes was used to survey the polymorphisms 
present at HvFTl, 89 cultivars, mainly European, and 20 inbred 
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lines derived from Spanish landraces (Igartua et al., 1998). They 
were classified into 60 winter and 49 spring types, according to 
their genetic constitution at vernalization and photoperiod genes 
(Table SI). These genotypes constitute a representative sample of 
barley germplasm available to European breeders. 

HvFTl phenotypic effect 

Several biparental populations used in past studies were 
reassessed to account for possible phenotypic effects of the poly- 
morphisms at HvFTl (Table 1). In some cases, these effects were 
already described in the references cited. In others, further geno- 
typing allowed a better resolution of the QTLs or the discovery 
of previously unknown polymorphism. QTL x QTL interaction 
analyses were done using the unbalanced analysis of variance 
option implemented in Genstat 14 (Payne et al, 2009), follow- 
ing a factorial model with the markers closest to the QTL peaks 
and "environment" as factors. The field experiments of the pop- 
ulation Beka x Mogador are explained in the publication by 
Cuesta-Marcos et al. (2008). 

Gene expression analysis 

Seven spring cultivars and two landrace-derived inbred lines, 
selected to represent HvFTl CNV types, were used for gene 
expression analysis. Also, eight doubled haploid (DH) lines of 
the population Beka x Mogador (Table 2) were used for this pur- 
pose. Several major flowering time genes were segregating in this 
population (Cuesta-Marcos et al., 2008). To focus only on varia- 
tion at HvFTl, DH lines were selected with spring (Beka) alleles 
at VrnHl, VrnH2, and PpdH2. Variation in HvCEN, the most 
important factor determining flowering time in this population, 
was also considered. We aimed at having two plants per HvCEN - 
HvFTl haplotype, but found only one for the Mogador-Mogador 
class (an extra plant was allocated to the Beka-Mogador class). 
PpdHl was not segregating in this cross. Both parents carry the 
recessive, long-day insensitive allele. 

SEQUENCE POLYMORPHISMS AT MAJOR GENES 

DNA sequence polymorphisms for the 109 accessions were 
screened with allele-specific primers of candidate genes. VrnHl 
was scored based on the size of the first intron of its candidate 
HvBM5A (Yan et al, 2003; von Zitzewitz et al, 2005). Alleles were 
classified according to Hemming et al. (2009); VrnH2 was evalu- 
ated as presence of HvZCCT-Ha and HvZCCT-Hb (Karsai et al., 
2005). PpdHl was genotyped using SNP22 in the CCT (Constans, 
Constans-like, TOC1 ) domain of its candidate gene HvPRR7, after 
digestion with BstU I (Turner et al., 2005). PpdH2 was scored as 
presence of the HvFT3 gene as reported by Casao et al. (2011). 
Regarding VrnH3, two indels in the promoter and two SNPs in 



the first intron of the HvFTl gene were assessed (Casas et al., 
2011). HvCEN, candidate gene for Mat-c or Eam6 was par- 
tially sequenced in 24 genotypes. The haplotypes are identified 
as reported by Comadran et al. (2012). 

It was not possible to assess HvFTl polymorphism directly at 
the Beka x Mogador DH population (Cuesta-Marcos et al., 2008) 
because the sequences of the parental alleles were conserved and 
they differed only in copy number (2 Beka, 1 Mogador). The 
population was reanalyzed based on the new information found 
(parents polymorphic at HvFTl for CNV). Two microsatellite 
markers, in the HvFTl region, were mapped (EBmac0603 and 
AF022725A), to provide better resolution of the flowering time 
QTL found in this region. 

HvFTl CNV ANALYSIS BY qPCR 

Genomic DNA was isolated from frozen barley leaves using the 
NucleoSpin Plant II kit (Macherey-Nagel, Germany) and used as 
template for CNV analysis by qPCR in an ABI 7500, essentially 
as described by Nitcher et al. (2013) with some modifications. 
Briefly, 100 ng genomic DNA were mixed with 2 |xM of each 
primer and 10 ul of 2X Power SYBR Green Supermix (Thermo 
Fisher Scientific, Waltham, MA). The PCR program comprised 
lOmin at 95°C, 40 cycles of 10 s 95°C and 50 s 60°C, and a 



Table 2 | Haplotypes for major flowering time genes and markers 
closest to QTLs (see Cuesta-Marcos et al., 2008) for selected doubled 
haploid (DH) lines of the population Beka x Mogador. 

DH VrnH2 VrnHl PpdH2 HvCEN HvFTl FLE + AE* 

line HvZCCT HvBM5 HvFT3 Bmac132 E41M47_e 



1840 B B B B B 42 

1892 B B B B B 41 

1927 B B B B M 43 

1933 B B B B M 44 

2011 B B B B M 44 

1873 B B B M M 30 36 

1837 B B B M B 31 37 

2009 B B B M B 40 46 



Alleles conferring lateness are highlighted in gray. Dates for developmental 
stages in the pot experiment used for gene expression analysis are also 
included. 

B, Beka; M, Mogador 

f FLE: days from sowing to 50% of plants with flag leaf expanded. 

*AE: days from sowing to 50% of plants with emerged awns. 

"more than 48 days, as the experiment was terminated at that date and the 

plants had not reached awn emergence yet. 



Table 1 | Biparental populations analyzed in this study. 



Population 

Henni x Meltan 
SBCC154x Beatrix 
Beka x Mogador 



Type Number of lines HvFTl polymorphism References 



DH 
DH 
DH 



118 
168 
120 



Promoter 
Intron and CNV 
CNV 



Borras-Gelonch et al., 2010 
Unpublished 

Cuesta-Marcos et al., 2008 



Present study 

Genotyping of HvFTl promoter 
QTL analysis, genotyping of HvFTl 
New markers in HvFTl region and CNV. 
QTL interaction 
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melting curve stage. Number of copies of the first exon of HvFTl 
was tested in all 109 genotypes. Morex was selected as the cali- 
brator genotype and SNF2 as the housekeeping gene (Yan et al., 
2002). Then, in a subset of lines, two other HvFTl regions 
(promoter and exon 3) and three other genes close to HvFTl 
(UCW118, UCW123, and UCW120) were tested as reported by 
Nitcher et al. (2013). Efficiency for each primer pair was obtained 
by serial dilutions of barley genomic DNA and it was taken into 
account for CNV calculation (Weaver et al, 2010). Efficiencies for 
SNF2, HvFTl -promoter, HvFTi-exonl, HvFTi-exon3, UCW118, 
UCW120, and UCW123 were 0.95, 0.96, 0.92, 0.95, 0.96, 1.03, and 
0.85, respectively. 

GROWTH CONDITIONS FOR EXPRESSION STUDIES 

Barley plants used for expression analysis were grown in pots of 
11.5 (diameter) x 16.0 (height) cm with a mix of peat, sand, and 
perlite. Six seeds of one genotype were sown per pot. After emer- 
gence, they were thinned to three seedlings per pot. The plants 
were grown in a growth chamber, under long-day conditions 
(16 h light, 250 u,E mT 2 s _1 , 20°C, 60% relative humidity/8 h 
dark, 16°C, 65% relative humidity) for 7 weeks. There were three 
pots per genotype, which were used as replicates. 

Two experiments were carried out with samplings at two dif- 
ferent times. In the first experiment, leaf tissue was harvested 
in the middle of the light period, after 8 h light, as reported by 
Kikuchi et al. (2009). In the second experiment, to maximize cir- 
cadian expression of HvFTl (Turner et al., 2005) harvesting took 
place 2h before dark, after 14 h light. In both cases, leaf tissue 
(last expanded or flag leaf) was harvested and frozen immedi- 
ately in liquid nitrogen before tissue homogenization (Mixer Mill 
model MM 400, Retsch, Germany). At each sampling time, three 
samples were analyzed per treatment and genotype. Each sample 
came from a different plant and pot. 

mRNA EXTRACTION. cDNA SYNTHESIS, AND GENE EXPRESSION 
ANALYSIS 

For qPCR analysis of HvFTl expression levels, 1 u.g of total RNA 
(purified using the NucleoSpin RNA Plant kit, Macherey-Nagel) 
was transcribed to cDNA by using the Superscript III reverse 
transcriptase and 2.5 u,M poly(dT)2o primer according to the 
manufacturer's instructions (Invitrogen). The reaction mixture 
for qPCR and the PCR program have been previously described. 
cDNA was quantified using a Nanodrop system (Thermo Fisher 
Scientific) and equal amounts were used for all samples. Actin was 
selected as the housekeeping gene (Trevaskis et al., 2006), and 
expression of HvFTl and HvCEN was analyzed using the same 
primers as in Yan et al. (2006) and Comadran et al. (2012), respec- 
tively. Efficiencies for Actin, HvFTl, and HvCEN, were 0.97, 0.86, 
and 0.96, respectively. 

RESULTS 

SURVEY OF HvFTl CNV 

Winter and spring genotypes were classified as follows: winter 
genotypes carry a functional VrnH2 allele and a winter or inter- 
mediate allele in VrnHl (wild type vrnHl, VrnHl-6, or VrnHl-4). 
Spring or facultative lines have been classified as those with a 
spring allele in VrnHl or lines in which VrnH2 is absent. With 



regard to HvFTl, there were accessions representing all possi- 
ble combinations of intron and promoter polymorphisms, both 
within the winter and spring groups (Table 3). The combination 
of the "late" promoter with the "early" intron, however, was the 
most frequent. HvFTl exon 1 copy number ranged from 0.35 ± 
0.13 in the winter cultivar Igri to 5.12 ± 0.28 in the spring cultivar 
Zaida (Figure 1). The results for several genotypes were inter- 
mediate between two classes, and were assigned to classes based 
on pedigree information when possible. Genotypes with several 
copies of HvFTl presented allele combinations typical of both 
winter and spring cultivars at all the Vrn and Ppd genes, with one 
exception, VrnH2. No genotype with the VrnH2 gene present had 
more than one HvFTl copy. 

It is remarkable that all of the winter barley cultivars, i.e., with 
a dominant VrnH2 allele, had only one copy of the first exon. 
Barberousse, the last winter genotype in Figure 1, had 1.50 ± 
0.01 copies, which we considered as a single copy. On the other 
hand, 28 out of the 49 spring cultivars analyzed contained more 
than 1.70 copies of exon 1 of HvFTl , which we have considered as 
multiple copies. 

The study was extended to analyze CNV in other areas within 
or near the HvFTl gene, as in Nitcher et al. (2013). Thus, 
as well as qPCR primers for HvFTl exon 1 (amplifying the 
region of +200 to +293 bp downstream from the ATG start 
codon), we used primers for UCW118 (nearest known gene in 
the flanking Morex BAC 455J22 upstream from HvFTl), HvFTl 
promoter (—727 to —656 bp upstream from the ATG), HvFTl 
exon 3 (+778 to +902 bp downstream from the ATG), UCW123 
(+6.6kbp downstream from HvFTl), and UCW120 (nearest 
known gene in the flanking Morex BAC 761F04 downstream 
from UCW123). This analysis was carried out in 10 barley vari- 
eties from Northern Europe (Figure 2A), 18 Spanish landraces 
(Figure 2B), and in another nine varieties from diverse origins 
(Figure 2C). Although all genotypes were not analyzed for all the 
genes, several results merit further attention. Two main patterns 
of HvFTl promo ter/exon 1/exon 3 copy number were found in 
Northern European barleys: Asplund, Olli, Herse, Stella, Tammi, 
and Maskin contained one single copy of the promoter and mul- 
tiple copies of exon 1 and exon 3, whereas Henni, Meltan, Pallas, 



Table 3 | Number of spring and winter accessions classified according 
to HvFTl haplotypes defined by polymorphisms at the promoter 
(indel 1- indel 2), first intron and number of copies of exon 1. 

Promoter Intron Number of copies 

1 >1 
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AG 
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135-146 
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139-142 


TC 
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7 
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FIGURE 1 | Determination of copy number variation for the first exon of HvFTl in winter (A) and spring or facultative (B) barleys. Bars represent 
means ± s.e.m. Morex was used as the calibrator genotype (copy number = 1). SBCC, Spanish Barley Core Collection. 



and Juli contained multiple copies of promoter, exon 1, and exon 
3. These different gene structures may affect gene functionality, 
as we will discuss later. CNV in other genes around HvFTl was 
only analyzed in 11 accessions. UCW118 (next to the promoter), 
and UCW123 and UCW120 (next to the exon 3) were similar to 
the corresponding CNV in the adjacent regions, with the excep- 
tion of Tammi, which contained multiple copies of UCW118 but 
only one promoter. We do not provide results for UCW120 in 
Meltan since the primers gave no amplification. Regarding the 
Spanish lines, only 4 out of 18 genotypes analyzed (SBCC157, 
SBCC154, SBCC133, and SBCC035) contained multiple copies of 
the promoter, exon 1 and exon 3. For the other genotypes ana- 
lyzed, only Beka and Dobla presented several gene copies in the 
HvFTl promoter, exon 1, and exon 3. 

EFFECT OF HvFTl POLYMORPHISM IN POPULATIONS 

We have reassessed several populations to illustrate the effect of 
sequence and CNV polymorphism at HvFTl . 

The population Beka x Mogador (Cuesta-Marcos et al., 2008) 
is a spring x winter population with very large flowering time 



variation. It has been reassessed because we have found that, 
although the sequence of HvFTl is conserved in both geno- 
types (1 distal mismatch in 2547 bp sequenced), it still presents 
CNV polymorphism: Beka has two copies of the gene, whereas 
Mogador has one. A small flowering time QTL in the vicinity of 
HvFTl was already reported, although the region was not well 
covered with markers in the original study. Close markers were 
now identified and genotyped in the population to increase cov- 
erage. The QTL already detected on 7HS in this cross was made 
more conspicuous with the markers introduced for this study 
(Figure 3), with the peak hinting at the HvFTl position, and a 
high significance [-loglO(P) above 11]. In this case, Beka con- 
tributed the early allele, about 2.2 days earlier than the Mogador 
allele. If this QTL is truly due to the effect of HvFTl, it must be 
caused by differential effect of the number of copies. The pos- 
sible effect of HvFTl CNV on the vernalization mechanism, as 
described by Nitcher et al. (2013), should have been evident in 
this population as an interaction of the QTL found at the VrnHl, 
VrnH2 regions with the QTL at the HvFTl region, in trials with- 
out enough vernalization. This could have occurred in the three 
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FIGURE 2 | Analysis of copy number variation for HvFTH promoter, exon 
1, and exon 3), and for regions upstream [UCW118) or downstream 
{UCW123 and UCW120) of HvFTl on barley genotypes of Northern 



European (A), Spanish (B), and other origins (C). Bars represent means ± 
s.e.m. Morex was used as the calibrator genotype (copy number = 1 ). SBCC, 
Spanish Barley Core Collection. 



late sowings, in which the number of cooling degree days (a mea- 
sure of vernalizing potential) was much lower than at the fall 
sowings (Table SI in Cuesta-Marcos et al., 2008). The six field 
trials were reanalyzed for interactions between the five major 
QTL {VrnHl, VrnH2, HvCEN, PpdH2, HvFTl), using the closest 
marker to each peak, and dividing the six trials into fall sow- 
ings (November) and winter-spring sowings (late February to late 
March). The effect of HvFTl was clear across all trials, but did not 
interact much with other genes. There was just one possible inter- 
action of HvFTl with VrnHl but in the fall-sown trials and not 
in the late sown trials. It was caused by a significant difference for 
the HvFTl alleles (with Beka, the one with two copies, inducing 
earliness) only in the presence of the winter (Mogador) allele at 
VrnHl (Table 4). The expected interaction between vernalization 
and HvFTl at the late-sown trials, however, was not detected at 



all. Therefore, we can conclude from this result that the possible 
promoting effect of the double HvFTl gene of parent Beka was 
not strong enough to affect the vernalization requirement of the 
winter lines of this population. 

SBCC154 x Beatrix is a cross of two spring genotypes. The pop- 
ulation was described by Hofmann et al. (2013), but the flowering 
time data have not been reported yet. SBCC 154 has four copies of 
HvFTl, whereas Beatrix has only one. Both have the putatively 
late (139-142) promoter (Nitcher et al., 2013), whereas Beatrix 
has the early intron (TC) and SBCC154 the late one (AG). The 
HvFTl marker detected a significant QTL with an effect of 2.5 
days, with Beatrix as the early allele (Figure 4). 

The population Henni x Meltan (Borras-Gelonch et al., 2010), 
a cross of two spring cultivars, was an example of polymorphism 
just at the promoter. Each genotype has four copies of HvFTl 
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FIGURE 3 | Multi-environment QTL analysis for days to flowering from 
January 1st for six field experiments carried out with 120 doubled haploid 
lines and the parents of the population Beka x Mogador. The peaks above 
the threshold (dashed line) indicate presence of QTL significantly affecting the 
trait. Data taken from Cuesta-Marcos et al. (2008), enriched with new markers 



on the 7HS chromosome. In the lower part of the figure, field trials are coded 
with Au (autumn sowing), Wi (winter sowing), or Sp (spring sowing) and two 
digits for the year; the colored dashes indicate the extent of the QTL and its 
direction: blue means that the early allele came from Mogador, yellow-brown 
from Beka, with intensity proportional to the size of the effect. 



and, although they may have just two copies of the promoter; 
the parents present the same number of copies across the whole 
gene. The only polymorphism found was at the promoter, with 
Meltan carrying the early promoter (135-146) and Henni the late 
one (139-142). Confirming this expectation, the early allele of the 
QTL in Borras-Gelonch et al. (2010) was contributed by Meltan, 
the difference being 23°C d (around 1-2 days). 

EXPRESSION ANALYSIS OF HvFTl IN SELECTED SPRING BARLEYS 

After finding CNV variation for HvFTl, the next step was to 
evaluate the effect of CNV variation on gene expression. If the 
effect of several copies of the gene was always as large as reported 
by Nitcher et al. (2013), then it should be detectable as a large 
increase of gene expression and a very early phenotype, over- 
riding the effect of any other polymorphisms at HvFTl. This 
hypothesis was tested by evaluating gene expression on a set of 
spring and facultative cultivars representative of different HvFTl 
copy number alleles: Morex (reference genotype for one copy of 
HvFTl) and Beatrix as single HvFTl copy number, SBCC154, 
SBCC157, Dobla, Beka, Pallas, and Juli as examples of multi- 
ple promoter and exon 1 copies, and Tammi as representative 
of the genotypes with one promoter but multiple exon 1 copies 
(Table SI and Figure 5A). HvFTl expression was found, in gen- 
eral, to increase according the developmental stage of the plants, 
although the expression levels differed widely among genotypes. 
At the second sampling date all genotypes, except Dobla, had not 
reached awn emergence yet. At this point, SBCC154, SBCC157, 
and Dobla, the earliest genotypes, displayed the highest HvFTl 
mRNA levels (awn emergence dates are included in Figure 5 A). 



Apart from this observation, we did not observe much corre- 
lation between HvFTl expression and time to awn appearance. 
Interestingly, these three varieties had the dominant allele for 
PpdHl, whereas the rest had the recessive allele at this gene. Also, 
we could not find a clear correlation between number of HvFTl 
copies or sequence polymorphism and HvFTl mRNA expression. 
For example, Juli and Tammi contained multiple copies of HvFTl 
(four promoters and four genes in Juli, one promoter and four 
genes in Tammi), but both of them showed lower HvFTl expres- 
sion than SBCC154, SBCC157, or Dobla, at the same sampling 
date. The influence of the difference in HvFTl promoter copies 
on HvFTl expression observed in Tammi with respect to Juli and 
Pallas will be discussed later. Thus, we can conclude that CNV has 
not a prevailing effect on HvFTl expression, over other types of 
polymorphisms at the same gene. It is not the main factor con- 
trolling HvFTl expression and flowering time (measured as awn 
appearance), and it depends largely on the genetic background, 
and polymorphisms at other genes, especially PpdHl, as has been 
reported previously in the literature in studies done before taking 
into account CNV (Turner et al., 2005; Hemming et al., 2008). 

EXPRESSION ANALYSIS OF HvFTl IN SELECTED LINES OF THE BEKA x 
MOGADOR POPULATION 

In a second gene expression experiment, eight DH lines of the 
Beka x Mogador population were analyzed, together with some 
genotypes in common with the first experiment (Figures 5B,C). 
This population was found to contain a QTL for flowering time in 
HvFTl gene, as shown above. Individuals were selected according 
to CNV in HvFTl: Beka contained two copies while Mogador 
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contained only one. qPCR analysis were performed as before, but 
this time the material was harvested 2 h before dark. Sampling 
time actually had a major effect on the detection of HvFTl 
expression: under these conditions, HvFTl mRNA levels were, in 



Table 4 | Probabilities of significance from three analyses of variance 
of flowering date recorded at six field trials, three fall-sown and three 
winter- or spring-sown in the Beka x Mogador population. 



Source 


All trials 


Fall-sown 


Winter- spring-sown 


CPROB* 


Trial 


0.000 


0.000 


0.000 


VrnH1 


0.000 


0 280 


0.000 


VrnH2 


0.000 


0.000 


0.000 


HvFT3 


0.000 


0.000 


0.000 


HvCEN 


0.000 


0.000 


0.000 


HvFTl 


0.000 


0.000 


0.003 


Trial. VrnHl 


0.000 


0.237 


0.000 


Trial. VrnH2 


0.000 


0.511 


0.000 


Trial. HvFT3 


0.092 


0.006 


0.632 


Trial. HvCEN 


0.000 


0.612 


0.001 


Trial. HvFTl 


0.766 


0.248 


0.698 


VrnH1.VrnH2 


0.000 


0.424 


0.000 


VrnH1.HvFT3 


0.489 


0.527 


0.644 


VrnH2.HvFT3 


0.947 


0.934 


0.811 


VrnHl HvCEN 


0.000 


0.008 


0.000 


VrnH2. HvCEN 


0.786 


0.072 


0.074 


HvFT3. HvCEN 


0.817 


0.305 


0.276 


VrnHl HvFTl 


0.097 


0.011 


0.699 


VrnH2. HvFTl 


0.564 


0.760 


0.652 


HvFT3. HvFTl 


0.330 


0.509 


0.477 


HvCEN. HvFTl 


0.109 


0.078 


0.548 



The sources of variance are the trials plus five markers close to QTL peaks rep- 
resenting five major flowering time genes: VrnHl (HvBM5), VrnH2 (HvZCCT), 
HvFT3 (Bmag382), HvCEN (Bmac132), and AFLP E41 M47_e (HvFTl). Threeway 
interactions are not shown, as none was significant. The analyses of variance 
were done on genotype means, taking as error the residual genotypic variance. 
'Conditional probability of significance for each term, when added to a full model 
with the rest of terms already included. P-values below 0.05 highlighted in bold 
type. 



general, clearly higher than at the first experiment (Figure 5A), in 
which samples were harvested 8 h after the commencement of the 
light period, instead of 14 h. The varieties used in the first exper- 
iment, which had a recessive ppdHl allele, were also included for 
comparison. Apart from the clear induction of expression from a 
dominant PpdHl allele, we were not able to establish a clear rela- 
tionship between HvFTl copy number, mRNA levels and flow- 
ering time for these plants (Figure 5B). HvFTl expression in the 
DHs was similar and stable during the experiment after 19 and 33 
days. No significant differences were observed in HvFTl mRNA 
levels between the two alleles as main effects, although there was a 
significant interaction (P = 0.022) between the alleles at HvCEN 
and HvFTl. At the second sampling date, the plants with the Beka 
allele in HvFTl showed significantly higher expression than plants 
with Mogador allele, only if the allele at HvCEN also came from 
Beka. The Mogador allele at HvCEN had a major effect on ear- 
liness, and these lines (1837, 1873, and 2009) flowered markedly 
earlier than the others (flag leaf unfolding dates are included in 
Figure 5C). 

Another interesting result of the second experiment was that 
Tammi showed the largest HvFTl expression after 33 days, and 
also flowered the earliest. It seems that Tammi has a more effec- 
tive HvFTl gene to promote flowering than the others. This could 
be related to the fact that Tammi, as indicated above, presented 
four copies of the HvFTl gene but only one copy of the pro- 
moter, whereas Juli and Pallas, for example, have four copies of 
both promoter and exon 1. 

EXPRESSION ANALYSIS OF HvCEN IN THE BEKA x MOGADOR 
POPULATION AND IN SELECTED SPRING BARLEY GENOTYPES 

As we have shown previously, CNV and different haplotypes 
for HvFTl are not sufficient to establish a clear relationship 
between them, HvFTl mRNA levels and awn appearance. For 
example, the presence of a dominant PpdHl allele was found 
to enhance HvFTl expression in different spring genotypes 
(Figure 5A). Additionally, CNV in HvFTl did not have a clear 
role to determine flag leaf unfolding in the Beka x Mogador 
DHs (Figure 5C). For these reasons, we decided to analyze the 
mRNA levels of another gene involved in the establishment 
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FIGURE 4 | QTL scan for flowering time at a field trial for a population of 168 doubled haploid lines from the cross SBCC154 x Beatrix. The figures 
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FIGURE 5 | HvFTl expression analysis by qPCR. In a first experiment 
(A), leaf tissue was harvested in the middle of the light period (after 
8h light) 14, 28, or 41 days after sowing. Days until awn appearance 
(triangles) are shown for the whole duration of the experiment (50 
days). HvCEN, HvFTl, and PpdH1 haplotypes for each plant are also 
indicated. In the case of HvFTl, polymorphisms for promoter (E, 
"early"; L, "late"), intron 1 (AG or TC) and copy number variation (for 



of flowering time, HvCEN, a paralog of FT and TFL1 in 
Arabidopsis (Kobayashi et al., 1999). We used the same cDNAs 
as for the previous expression studies (Figure 6). We included 
haplotype information for some genotypes about HvCEN in 
Table SI for comparison. As observed for HvFTl, our results 



exon 1) are shown. In the second experiment (B,C), harvesting took 
place 2h before dark (after 14 h light), 19 or 33 days after sowing. 
Days until full unfolding of the flag leaf (triangles) are shown for the 
whole duration of the experiment (50 days). As previously, HvCEN and 
HvFTl haplotypes for each plant are included. (C) Expression in 
selected doubled haploid lines of the Beka x Mogador population (BxM, 
see Table 2). Bars represent means ± s.e.m. 



suggest that HvCEN mRNA levels are not the main respon- 
sible to explain awn appearance (Figure 6). HvCEN expres- 
sion dramatically increased when samples were harvested 14 h 
after dawn instead of 8h (Figures 6 A-C), as we observed 
for HvFTl. Another remarkable observation is that HvCEN 
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FIGURE 6 | HvCEN expression analysis by qPCR. HvCEN mRNA levels 
were quantified using the same conditions as for Figure 5. (A) leaf tissue 
was harvested in the middle of the light period (after 8 h light) 14, 28, or 41 



days after sowing. (B,C), harvesting took place after 14h light, 19 or 33 days 
after sowing. (C) Expression in selected doubled haploid lines of the Beka x 
Mogador population (BxM, see Table 2). Bars represent means ± s.e.m. 



expression also changed during the development, in disagree- 
ment with previously observed results for HvCEN expression 
(Comadran et al., 2012). For example, Morex contained signif- 
icantly higher HvCEN mRNA levels than Tammi, but Morex 
did not flower after 50 days and awns in Tammi appeared after 
33 days. 



DISCUSSION 

DISTRIBUTION AND PHEN0TYPIC EFFECT OF HvFTl POLYMORPHISMS 
OVER BARLEY GERMPLASM 

In this paper, we analyzed the extent of HvFTl CNV and its effect 
on HvFTl expression in more than 100 different spring and win- 
ter barley genotypes, mainly European, including some landraces 



Frontiers in Plant Science | Plant Genetics and Genomics 



June 2014 | Volume 5 | Article 251 | 10 



Loscos et al. 



HvFTl CNV and expression 



from the Spanish Barley Core Collection (SBCC) (Igartua et al., 
1998). Recently, Nitcher et al. (2013) demonstrated that the 
FT1 allele present in the barley genetic stock BGS213, which 
carried several copies of HvFTl, showed earlier transcriptional 
up-regulation of FT1, and was associated with a dominant spring 
growth habit. Cultivar Tammi, from Finland, was the donor of the 
mutation present in that genetic stock (Franckowiak and Konishi, 
1997). This allele is not frequent and was reportedly found only 
in spring cultivars originating from regions of extremely high lati- 
tude or high altitude (Takahashi and Yasuda, 1971). We confirmed 
the finding by Nitcher et al. (2013) that the Tammi/BGS213 allele 
was characterized by having a single copy of the promoter and 
several copies of the transcribed region. This kind of allele was 
present in a group of 6 six row Northern European cultivars, some 
of them related by pedigree (Baumer and Cais, 2000; Chiapparino 
et al, 2006). The origin of this allele could be traced back to prim- 
itive cultivar Asplund, which is in the parentage of Tammi and 
Herse, both showing alleles with multiple exon copies and one 
promoter (Figure 2). This allele seems to be particularly benefi- 
cial at high latitudes because both parents of Tammi (Asplund 
and Olli) come from the Northern limits of the barley cultivation 
range (Manninen and Nissila, 1997). 

We found several copies of HvFTl in accessions from appar- 
ently disconnected germplasm groups, like Scandinavian and 
Spanish landraces. Although the sample examined is not enough 
to derive definitive conclusions, the distribution of CNV alleles 
over type of growth habit haplotypes is intriguing. Multiple copies 
at HvFTl occurs only in genotypes that do not have winter growth 
habit. This suggests a possible disadvantageous agronomic effect 
of the presence of several copies of the gene in interaction with 
VrnH2. A possible mechanism that provides a plausible expla- 
nation for this is that the increased copy number of HvFTl is 
epistatic to winter alleles for VrnHl and VrnH2 (as reported by 
Nitcher et al., 2013). Munoz-Amatriain et al. (2013) have recently 
shown that high levels of CNVs are found in the barley genome, 
around 9.5% in coding genes, at increasing frequencies as distance 
to the centromere increases. This is a widespread phenomenon, 
contributing to phenotypic variation in barley. HvFTl CNV may 
be neutral in spring genotypes, but breeders could have selected 
against this duplication when breeding winter barley to avoid 
early transition to reproductive growth and exposure of repro- 
ductive tissues to damaging low temperatures. However, we have 
shown that this epistatic mechanism is not universal, as it does 
not occur in the Beka x Mogador population. 

The reassessment of three biparental populations representing 
polymorphisms of different kind identified QTL in the region of 
HvFTl. Previous studies also shed some light on the possible phe- 
notypic effect of HvFTl. The QTL found in the cross SBCC145 x 
Beatrix (Ponce-Molina et al, 2012), whose peak coincides with 
the gene, must be due to polymorphisms at the promoter region, 
as both parents have one copy of the gene and same first intron 
sequence. The early allele was contributed by SBCC145, a Spanish 
landrace from the Canary Islands, which carries the "early" pro- 
moter. It must be pointed out, however, that SBCC 145 has a 
distinct promoter, with additional polymorphisms compared to 
others, similar to cultivars Dairokkaku (Casas et al., 2011) and 
Meltan. The cross SBCC 154 x Beatrix offered the opportunity to 



contrast two different polymorphisms (promoter for SBCC145 x 
Beatrix, first intron, and CNV for SBCC 154 x Beatrix) against a 
common parent (Beatrix). This population was tested at a sin- 
gle trial sown in autumn, in the same field as the population 
SBCC145 x Beatrix, and the result was the presence of a QTL, 
exactly at the marker used to genotype HvFTl, with SBCC154 
as the late allele. In this case, high copy number (SBCC154) was 
late compared to one copy, but this could be due to an effect 
of SBCC154 carrying the "late" intron. Therefore, any acceler- 
ation of flowering that might be produced as a consequence of 
high copy number was secondary to the lateness associated to the 
sequence polymorphism. It must be noted that PpdHl was segre- 
gating in this population. However, flowering occurred too early 
to allow for a significant effect of PpdHl , which was not detected 
as QTL. 

In the cross SBCC016 x Esterel (Casas et al., 2011), each geno- 
type had one copy of HvFTl, Esterel has the "early" intron (TC) 
and SBCC016 the "late" one (AG). Accordingly, the early allele 
of the QTL, which also peaked at the gene itself, was contributed 
by Esterel. Hayes et al. (1993) and Borras-Gelonch et al. (2012) 
reported a QTL in this region for the population Steptoe x Morex. 
Although the size of the effect was not large, in both studies the 
Morex allele was significantly later. Both parents have just one 
copy of the gene, and a mixture of late and early polymorphisms 
at the promoter and the intron. It seems that the intron effect is 
stronger, as Steptoe carries the early (TC) allele at this position. 

A QTL in the region of HvFTl was also detected in the classi- 
cal studies carried out in the Igri x Triumph population (Laurie 
et al., 1995), with the Triumph allele conferring earlier flower- 
ing, although the nucleotide sequences were identical (Yan et al., 
2006). CNV has been recently identified in this population, with 
two copies of the HvFTl gene in Triumph (R. Nitcher, personal 
communication). This result would match our findings in the 
Beka x Mogador population. Nevertheless, we could not repro- 
duce this result, since the Triumph seed we analyzed had only one 
copy of the gene. We cannot discard that the samples analyzed 
are different, because this cultivar is known under two different 
names, Trumpf in Germany and Triumph in the UK, where it 
was reselected from somewhat heterogeneous seed (van Harten, 
1998). 

In summary, these findings reported in the literature, com- 
bined with the results presented in this study, reveal that detection 
of flowering time QTL in the region of HvFTl in biparental 
populations representing all kinds of polymorphisms at HvFTl 
(promoter, first intron, and CNV) is common. There is no func- 
tional proof that this gene is responsible for all the QTL, but it is a 
good candidate. An alternative explanation could be the presence, 
at least in some cases, of an additional flowering time gene closely 
linked to HvFTl. In any case, we have to wait until there are 
either functional proofs or increasing evidence from other pop- 
ulations to declare that HvFTl is the responsible for the variation 
observed. We expect that the catalog of polymorphisms presented 
in this study will help other researchers to contribute information 
to clarify the issue. 

The survey of HvFTl polymorphisms allows us to conclude 
that HvFTl (VrnH3) is far from being effectively fixed in culti- 
vated barley, as stated up to now (Stracke et al., 2009; Comadran 
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et al, 2012). This statement is probably true if one only considers 
the allele responsible for the huge phenotypic effect observed by 
Nitcher et al. (2013) on a winter barley, which confirmed obser- 
vations of the seminal work by Takahashi and Yasuda (1971). The 
effects that we have found in different populations (this study), 
association studies (Casas et al., 201 1) and also reported in the lit- 
erature (previous paragraphs) point at smaller phenotypic effects 
that are linked to all types of polymorphisms in this gene (pro- 
moter, first intron, CNV). These effects, however, cannot be easily 
extrapolated to different genetic backgrounds. 

CNV AND GENE EXPRESSION 

The commencement of the reproductive stage in barley and the 
duration of the time period until flowering are controlled by a 
variety of factors that act interactively. In temperate cereals like 
wheat and barley, flowering is promoted by long days. In barley, 
the up-regulation of VRN3/FT1 appears to be the main trig- 
ger for the initiation of flowering (Faure et al, 2007), although 
an alternative view is that its role may be to accelerate inflores- 
cence development and reduce the time taken from double-ridge 
to head emergence (Sasani et al., 2009). Its expression is tightly 
regulated, repressed by VRN2 under long days (Hemming et al., 
2008), which actually competes with CO to bind to protein com- 
plexes that activate -FT (Li et al., 2011). Also in response to long 
days, FT1 expression is promoted by CO-like proteins (Li et al., 
2011; Campoli et al., 2012) and the pseudo response regulator 
PpdHl (Turner et al., 2005; Faure et al, 2007). This promotion 
may occur by interaction of CO with the promoter of FT, as hap- 
pens in Arabidopsis (Tiwari et al., 2010; Li et al., 201 1; Andres and 
Coupland, 2012). 

The findings of this study add to mounting evidence pointing 
at a complex control of the timing of head emergence that cannot 
be easily reduced to a simple scheme. On one hand, HvFTl dis- 
plays a variety of polymorphisms in regions that are compatible 
with a regulatory role. On the other hand, its expression pattern is 
compounded by its nodal position, downstream of the vernaliza- 
tion and long day pathways, whose genes also have large influence 
on duration of plant development. 

The first gene expression experiment results confirmed the 
induction of HvFTl expression by the long day pathway, irrespec- 
tive of the polymorphisms present at HvFTl. The genotypes that 
reached first awn appearance in that experiment were all those 
having a sensitive/dominant PpdHl allele, particularly at the sec- 
ond sampling date. This date is probably the most meaningful 
because all PpdHl genotypes reached awn appearance just a few 
days later. This result agrees with other reports where HvFTl 
expression in ppdHl background was markedly lower (Turner 
et al, 2005; Hemming et al, 2008; Kikuchi et al, 2009). Therefore, 
we can conclude that gene duplication of HvFTl does not always 
ensure higher expression and, in any case, the scale of its effect 
is minor compared to the induction by PpdHl . The case of high 
HvFTl expression in presence of ppdHl described in Nitcher et al. 
(2013) seems an exception, probably due to a unique genetic 
background. 

Our results confirm that the triggering of events at the meris- 
tem leading to flowering is not determined just by HvFTl expres- 
sion. There must be additional factors that probably need the 



presence of HvFTl product to interact with. We cannot rule out, 
however, that copy number is related to a dosage effect, precisely 
by the different genetic background in each genotype, which 
may lead to differences in the induction of HvFTl. A benefi- 
cial increase in dosage is actually one of the evolutionary forces 
explaining the conservation of gene duplications, particularly for 
genes that mediate the interaction between the organism and the 
environment, as reported for Ppd-Bl and Vm-Al in wheat by 
Diaz et al. (2012), or for genes with dosage-sensitive functions, 
owing to protein-protein interaction (Innan and Kondrashov, 
2010). HvFTl corresponds to the first class, and may as well be 
part of the second. CNV is widespread in plants (Zmienko et al., 
2014), and certainly in barley (Munoz-Amatriain et al., 2013), but 
its effects vary in each case. For instance, in barley, the powdery 
mildew resistance allele mlo-11 acts by disrupting its expres- 
sion through accumulation of non-functional copies of the gene 
upstream of the wild-type copy (Piffanelli et al, 2004). On the 
other hand, the effect of the freezing tolerance locus FrH2 in bar- 
ley depends on the number of CBF genes transcripts produced, 
that are proportional to the number of copies present in the gene 
cluster identified as responsible for this QTL (Pasquariello et al., 
2014). The plant immunity locus GER4, also in barley, is also a 
cluster of tandemly duplicated genes. In this case, the enhanced 
transcript dosage was proposed as the evolutionary driving force 
for the local expansion and functional redundancy of the locus 
(Himmelbach et al, 2010). 

The results of Nitcher et al. (2013) and the earliness induced 
by the double copy "Beka" allele in the Beka x Mogador popula- 
tion point in the direction of a dosage effect. However, the direct 
comparison of HvFTl expression figures for different alleles is 
hampered by the large influence of the alleles at PpdHl on HvFTl 
and at HvCEN on flowering time. The analysis of the eight sister 
lines of the Beka x Mogador population indicates a slightly ear- 
lier flowering and higher HvFTl expression of lines with the Beka 
allele only in presence of the Beka allele at HvCEN (spring, Hap 
III), but the evidence is too small to declare that we have found a 
clear dosage effect. 

These results led us to speculate that the most efficient ver- 
sion of the barley gene could be one with a single conserved 
promoter, and variable number of transcribed regions, as seen 
in cultivar Tammi, or one promoter and one copy of the gene. 
Other versions of the gene, with several promoters and tran- 
scribed regions may take longer to be induced. Indeed, we do 
not know whether the multiple copies of HvFTl are functional 
or not. They seem to be expressed in Tammi/BGS213, but there 
is no evidence for other alleles. If they are, the presence of sev- 
eral full and functional copies of the gene may affect differently 
the expression of the gene, depending on the concentration of 
the promoting signal. If it is low, its dilution over several copies 
may delay transcription, whereas an abundant promoting signal 
would enhance transcription proportionally in alleles with sev- 
eral copies. The special case of Tammi/BGS213 indicates either 
that the transcription of several copies is triggered by a single 
promoter, or that the additional copies of the gene are placed 
under the control of other promoters and thus escape the regula- 
tory control of the wild-type gene. For other multiple copy alleles, 
with equal number of promoters and transcribed regions, we can 
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speculate that the copies are full and perfect replicates of the orig- 
inal gene, and that the expression of them all is affected by the 
same processes. 

The difference in overall gene expression between the two 
experiments agrees with the reports on pattern of diurnal oscil- 
lation of FT1 expression during the day (Turner et al., 2005; 
Kikuchi et al, 2009). Actually, the same observation can be 
made for HvCEN (Figure 6). This information, combined with 
the sequence analogy between the two genes, suggests also that 
HvCEN could be subject to circadian rhythm. 

INTERACTIONS OF HvFTl 

The DH lines (Figures 5C, 6C) presented very different dates 
to reach a developmental stage (flag leaf unfolding in this case) 
in the presence of rather similar levels of HvFTl expression. 
As pointed out before, in this case the cause for differences 
in development seems to lie on the allele present at HvCEN. 
It is clear that accelerated development in the presence of the 
Mogador HvCEN allele (Hap II) did not depend on higher 
expression of HvFTl. Therefore, early development in these 
lines must be caused by some mechanism acting in parallel to 
HvFTl, either through protein-protein interaction with HvFTl 
or by a combined effect with HvFTl by focusing on the same 
targets. 

Actually, the interaction between HvFTl and HvCEN is by 
no means unexpected. HvCEN and HvFTl are probably paralogs 
as suggested by their high sequence identity (59%), and indeed 
have been annotated as members of closely related protein fam- 
ilies before (Higgins et al., 2010; Andres and Coupland, 2012; 
Comadran et al., 2012). Therefore, it is likely that their prod- 
ucts display similar interactive patterns with other proteins. There 
are many evidences in the literature pointing at the interaction 
of FT proteins with other proteins in the flowering promotion 
pathway. Ahn et al. (2006) reported that the family of small pro- 
teins coded by FT and homologous genes act as "either scaffolds 
or regulators of signaling complexes" in Arabidopsis. The bZIP 
transcription factor FD also seems to play a central role in flow- 
ering (Wigge et al., 2005). More recently, Jaeger et al. (2013) 
put forward a model to explain flowering control in Arabidopsis 
in which TFL1 (the product of another paralog gene) com- 
peted with FT to form the complex with the FD gene product, 
needed to trigger floral meristem specification. Some evidence for 
this kind of process in cereals was found by Li and Dubcovsky 
(2008), who detected interaction between the proteins of TaFT 
(the FT homolog of bread wheat) and FD-like2, and also in 
rice (Taoka et al., 2011), although in this case the interaction 
needed the intermediation of a third protein. In barley, how- 
ever, FD orthologs can only be predicted in terms of sequence 
similarity. 

A search in the Protein Data Bank revealed that, besides the 
FT-FD interaction, Arabidopsis homodimers of FT and also TFL1 
have been reported. Thus, it is also conceivable that proteins 
coming from genes as similar as HvFTl and HvCEN may form 
heterodimers as well, making the speculation about the interac- 
tion between the products of these two genes more plausible. The 
HvCEN polymorphism, which translates to a non-synonymous 
mutation (Prol35Ala), is located in a solvent-exposed protein 



loop (Comadran et al., 2012). However, with the structural 
evidence at hand, this loop does not directly participate in any 
dimeric interface and therefore nothing can be concluded about 
its molecular role in the interface. 

The role of HvFTl may actually go beyond the duration 
of developmental phases. There is recent evidence about dra- 
matic agronomic effects of the ortholog of FT in tomato, SFT 
in interaction with SP (SELF PRUNING, itself an ortholog of 
HvCEN and HvTFLl). SFT is in a "dose-dependent epistatic 
interaction" with SP which results in a modification of plant 
architecture that can be optimized to produce higher yields 
(Jiang et al., 2013). This dose-dependent action of SFT 
also opens the ground for speculation about possible dose- 
dependent action of multiple copies of HvFTl and its agronomic 
outcome. 

CONCLUSIONS AND FURTHER WORK 

The main conclusions of this study are: 

• A wide survey of barley germplasm revealed that HvFTl dupli- 
cation was only observed in spring and facultative barleys that 
do not possess a functional VrnH2 allele. 

• Two models of HvFTl duplication were observed, one that 
includes the promoter and the gene, the other only the tran- 
scribed region. Higher gene expression seems associated only 
to the second one. 

• There are flowering time QTL on the region of HvFTl in dif- 
ferent populations, representing all types of polymorphism at 
HvFTl , promoter, first intron and CNV. 

• Analysis of HvFTl expression and phenotypic effects showed 
that they depend on gene polymorphisms but also on genetic 
background. 

Plant breeders must be able to fully harness the development of 
cereal plants to be able to respond to the challenges of climate 
change. In this study, we present the case to state that HvFTl has 
been a hot spot to fine tune barley adaptation to environmen- 
tal conditions, and will have to be given due consideration by 
breeders to create future cultivars. 

The role of other genes in triggering flowering initiation, pos- 
sibly in interaction with HvFTl, seems a very rich area which 
deserves more research. In particular, further avenues to attempt 
to resolve the precise sequence or structural variation in HvFTl 
which causes early or late flowering,- particularly for the Tammi 
allele will require crossing different HvFTl alleles into a common 
genetic background, maintaining fixed PpdHl and HvCEN. This 
could be combined with BAC sequencing for the entire HvFTl 
region, for key alleles. 
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